Interchain coherence of coupled Luttinger liquids at all orders in perturbation theory 
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We analyze the problem of Luttinger liquids coupled via a single-particle hopping t± and introduce 
a systematic diagrammatic expansion in powers of t±. An analysis of the scaling of the diagrams 
at each order allows us to determine the power-law behavior versus t± of the interchain hopping 
and of the Fermi surface warp. In particular, for strong interactions, we find that the exponents 
are dominated by higher-order diagrams producing an enhanced coherence and a failure of linear- 
response theory. Our results are valid at any finite order in t± for the self-energy. 
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The problem of crossover from one to higher dimen- 
sions has recently received particular interest The 
puzzling question is whether and how an infinite row of 
one-dimensional chains of Luttinger liquids (LL) develops 
interchain coherence and possibly goes over into a Fermi- 
liquid state for an arbitrarily weak single-particle hop- 
ping t± <C Ep between the chains {Ep being the Fermi 
energy), and how this is affected by the correlation ex- 
ponent a flo[| of the unperturbed LL which parametrizes 
the intrachain interaction. No clear agreement has been 
reached so far and different analyses have provided con- 
flicting answers: (i) The system might go over to a Fermi- 
liquid state for arbitrarily small t±, or (ii) there might 
be a finite t± (depending on a) below which the sys- 
tems remains in a LL state and no coherence is devel- 
oped between the liquids 0, or finally (iii) one could 
have a different state possibly with gaps and long-range 
order. This ambiguity is essentially due to the fact that 
the perturbation introduced by t ± is relevant ^ || . This 
means that no matter how small ij_ , perturbative results 
in t±_ will be always inaccurate for sufficiently low values 
of some characteristic energy £ < £ = (t±/ Ep) 1 ^ 1 ^ ^. 
Here, £ is the largest energy scale between (i) the ex- 
ternal frequency w, (ii) the momentum measured from 
the Fermi surface \k — kp\ JhJ , and (iii) the temperature 
T. Indeed, in the interesting region £ ~ £q, where in- 
terchain coherence starts to set on [0, all perturbation 
terms are of the same order and no definitive predic- 
tion can be made about the nature of the ground state. 
An alternative approach to this topic is to consider the 
limit of strong forward scattering where one can use the 
higher-dimensional bosonization method or the Ward 
identities B. In this limit, one can show that a system 
always becomes a Fermi liquid in dimensions greater than 
one if the interaction is not too singular. 

In this Letter, we analyze how interchain coherence 
develops when t± is switched on and its a dependence. 
The striking point of the present work is that the results 
are valid at any finite order in t±_ for the single-particle 
self-energy, in contrast to previous work ||,||||||. This 
is important because it fixes some essentially exact re- 
sults to which approximate theories should be compared, 
although it is clear that an answer to the more inter- 
esting question "Fermi liquid or not" to this degree of 



accuracy is extremely difficult and will not be addressed 
here. In addition, we show that restricting to lowest- 
order terms may give incorrect results even in the small 
t± limit for a larger than a certain value. In particu- 
lar, when considering higher-order diagrams, interchain 
coherence is increased, and linear-response theory does 
not work. Specifically, we evaluate the exponent of the 
power-law behavior in tj_ (i) of the interchain hopping 
of a particle at the Fermi momentum, and (ii) of the 
Fermi surface (FS) warp. To calculate these quantities 
we extend a method developed in Ref. Jll]] to the case 
of coupled LL and rewrite diagrammatically the pertur- 
bation expansion in t± introduced in Ref. ||. We then 
consider the low-energy scaling behavior of the diagrams 
at all orders in ij_ and notice that the scaling is no longer 
homogeneous for a > a 2p (cf. also Refs. ||U). Finally, 
we use this scaling behavior to calculate the exponents 
mentioned above. Recent exact-diagonalization calcula- 



tions 1 12 1 are in excellent agreement with our predictions. 

We consider a system of N coplanar coupled identi- 
cal chains lying parallel to the x-axis with Hamiltonian 

H = T,i H LL(i) + Ei,V<r*-Uj J dx ^lA x ^)^a(x;j), 

where Hll(i) describes a (uncoupled) Luttinger liquid 
in chain i. t±ij is a hopping term between chain i and 
j, and tp rt(T (x; i) (V4 <r( x \ *)) i s the destruction (creation) 
operator for a right- (r = +1) or left-moving (r = — 1) 
fcrmion at site x in the chain i with spin a (we will also 
consider the spinless case) . The diagrammatic expansion 
in t around the atomic limit of the Hubbard model intro- 



duced by W. Metzner in Ref. 1 1 can be readily extended 
to the present case by considering perturbations in ij_ 
about the exactly-solvable Hamiltonian for the isolated 
Luttinger liquids ^iHll^i). We consider the diagram- 
matic expansion for the Green's function Q{_ku,iu; k±) 
expressed in the imaginary-frequency (iuj) and momen- 
tum representation, where k\\ and k± label the Fourier 
transforms of the x and of the i components, respectively. 
Unless otherwise specified, fcy will also implicitly con- 
tain the other in-chain quantum numbers r and a. The 
expansion of Q(ku,iu};kx) in ij_ is obtained by drawing 
all connected diagrams composed of two external lines, 
and of an arbitrary number of vertices connected with 
directed lines, each vertex having a number of entering 
lines equal to the number of leaving lines (see Fig. [l]). 
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One then labels each internal line (I) by a momentum 
(fcii^, fcx ) and a frequency variable (iu>^), and the 
external lines by fcii, k± and iuj. Each internal line (I) 
yields a factor ix[&x ] (the Fourier transform of t±y). 
Each vertex with 2n legs (2n > 2) acquires a contribution 
CjJOfen^W 1 ), • • • /^"W^^W 1 )', • • • k^'iuj^ 1 ), 
where C° is the n-particle cumulant of the isolated chain 
(for example, Cf = G°, the LL Green's function), and 
kn®iw( l \ kvS l ''iu>^' are the frequency and momenta 
(which are conserved at the vertex) of the entering and 
leaving lines, respectively (cf. |jfl]|5)). The introduction 
of higher-order cumulants and vertices is necessary, since 
the correlation functions of the t± = Hamiltonian do 
not satisfy Wick's theorem. Eventually, one has to mul- 
tiply each diagram by the usual symmetry and fermion 
factors and integrate (or sum) over all internal frequen- 
cies and momenta. This is a diagrammatic representation 
of the expansion introduced by D. Boies and coworkers 
in Ref. J5|. Similarly to conventional diagrammatic the- 
ory (cf. also Ref. p[ ), one can introduce a "self-energy" 
r(fc|| , iu>; k±) . The "inverse-self-energy" T(k\\,iw,kx) 
is then obtained by summing all one-particle-irreducible 
contributions to G(k\\ , iuj; k±) (see Fig. |l|c). One then 
obtains an analogue of the Dyson equation (cf. Fig. |l|b): 

g{k h iuj;k x ) = (r^n.iwifcx)- 1 -t±[k x ])~ 1 . 

At this point, it is clear that the approximation used 
in Refs. ( anc ^ others) consists in considering only 

the contribution to T from the first diagram (the dot a) 
in Fig. |l|(c) ["single-dot" approximation (SDA)] and thus 
using r = Q°. Furthermore, one can obtain the limit of 
infinite coordination D — > oo by summing an infinite se- 
ries of "bow" diagrams with "full" self-consistent Green's 
functions as internal lines (cf. Ref. [pT|). Nevertheless, 
we briefly review the results obtained by others [§,[|]8| 
within this SDA. Here, the system (both spinless and 
spinful) behaves as a Fermi liquid for any values of the 
correlation exponent a < 1 fiof , in the sense that there 
is a warped (i. e., fcx-dependent) FS fcy^fcjj, on which 
the Green's function has a (real) pole with finite weight 
at bj = and there is a finite region in momentum space 
around fc|| F [fcj_] where the pole survives, remains real, 
and continuously shifts to higher binding energies Jl3| . 
The quasiparticlcs are thus well defined everywhere on 
the FS (except at the special points for which tx[&x] = 0). 
Other poles Q do not describe true quasiparticles, since 
their imaginary part is always of the same order as their 
energy. 

However, the trouble starts to arise when one tries 
to consider higher-order diagrams. A simple power- 
counting argument indeed shows that a diagram of or- 
der tx™" 1 for T diverges like in the \aw-S 
limit. If one tries to keep the effective expansion pa- 
rameter t±Ep°'£ a ~ 1 small, then the tx[&x] term in the 
denominator of Q(k , iuj; k±) becomes much smaller than 
r(feii , iu>; Ax) -1 and the system behaves as a sum of un- 
coupled LL. On the other hand, as soon as one tries to 
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FIG. 1. Diagrammatic expansion in t± of the sin- 
gle-particle Green's function Q (gray box). Directed lines 
give a contribution £j_[&x]> and a dot with n entering and 
n leaving lines contributes a factor C° (n-particle cumulant). 
(a) Example of single-particle irreducible and reducible con- 
tributions to Q. (b) Dyson's equation for Q in terms of the 
inverse-self-energy V (gray disk), (c) Example of diagrams 
contributing to V (some of them are discussed in the text). 



approach the pole, the approximation fails and higher- 
order diagrams become of the same order of magnitude. 
For this reason, nothing reliable can be said about pos- 
sible Fermi-liquid behavior at this order of perturbation 
and it is mandatory to consider higher-order diagrams. 
In addition, the above power-counting argument does not 
hold for all values of a. For a greater than a certain 
ui v some diagrams turn out to diverge stronger at low 
frequencies (cf. also Refs. |],||||). For example, the con- 
tribution to t±T(xi — x 2 ;/cx) H] from the diagram 7 
(Fig. |l|c), with internal lines 3,4,5 taking the r indices 
+,—,—, respectively, is proportional to p5[ 
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(1) 



where |x| = x 2 + (|r| + a) 2 , argx = arg(x + it), and 
the exponent B = (1/K p — K p )/(2S). A naive dimen- 
sional analysis yields a term cx (t±/Ep) 4 |xi — X2| 2 4q , 
giving a contribution proportional to (t±/(Ep£ 1 ^ a )) 
to the Fourier transform t±T(kn, iui) k±). This is cor- 
rect for £ -c Ep if the integral does not diverge at 
small distances for a — > 0, i.e. as long as B < 1. 
For B > I [[l0| one picks up an a-dependent contri- 
bution a 2 ~ 2B ~ E F B ~ 2 , which must be balanced by 
£ 2 ~ 2B to give the correct energy dimensions. This pro- 
duces a contribution to t±T(k\\,iuj; k±) proportional to 

(t±/(E F £ 1_Q )) 4 (E F /£) 2B ~ 2 , i.e. a stronger divergence. 
A similar analysis shows that the 7-legs diagram 6 in 
Fig. ^(c) and its generalizations with 2n + 1 internal 
legs (n odd integer) produce a contribution to t±T pro- 
portional to (t±/(E%£ 1 - a )) 2n+2 (E P /e) %n ( B - 1 '>. In 
fact, this is the strongest low-energy divergence one 
can attain for a contribution to t±T at a given order 
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t± 2n+2 . The regime B > 1 corresponds to a > ct2 P with 
ct2p = (vS^ + l — i.e., where two-particle processes 
become more relevant than one- particle processes ■ 

We now use the above results to calculate the power- 
law behavior for small tj_ of the quantity [Q 



do; [0(fcjr, ioj; n) — G(kp,iuj] 0)] 



(2) 



Ank F is also equal to the expectation value of the single- 
particle hopping operator at kp in the two-chains case 
and is thus a measure of the coherence || of single- 
particle hopping. For k — kp and T = the only remain- 
ing energy scale that can be associated with £ is ui. For 
a < a.2 V all diagrams scale in the same way as discussed 
above. Therefore, at any order in the t± expansion and 
for uj <C Ep the difference between the Green's function 
of the two bands Q(kp 1 iuj; ^) — G(kp, iuj; 0) can be written 
(say, for Im iw > 0) as E^^cj^g [t^E F a {iu) a - 1 ], 
where g[x] is a scaling function. The change of variables 
lu = x ij_V(i-«) gives An kp oc (t^/ E F ) a / {1 ~ a ^ times 
a dimensionless integral, which converges at large x for 
a < 1/2. To assure convergence at low frequencies it is 
essential to expand the self-energy (r _1 ) in t^and not di- 
rectly the Green's function, in order to avoid an unphysi- 
caZlow-frequcncy (or low-T) divergence. For a > 1/2 the 
integral diverges for large x and one needs to introduce 
a large-energy cutoff yielding a dominant linear behavior 
Ank F oc t±/Ep. The same exponents would have been 
obtained by simply using the SDA for T. Notice also that 
this result is consistent with the one in Ref. Q. 

In the a > ct2 P regime, one cannot straightforwardly 
extend the above discussion since the diagrams give non- 
homogeneous contributions. If one stops the expansion in 
t± for the self-energy at an arbitrary finite order t± mo+1 
one obtains (restricting to dominant terms at each order) 



g(k F ,iuj;TT) 

mo 

= tx } J a r , 

m— — 1 



(3) 



((^) tt - 1 V*i)' 



iuj \ 
Ef~) 



(1-B)[m] 4 



where [m]^ is defined as 4 int((m — 2)/4) + 2, int be- 
ing the (lower) integer part, and the a m are constants 
(do = 1) possibly dependent on the spin and charge ve- 
locities. Q(kp, ioj; 0) _1 is obtained from the same expres- 
sion by replacing t± — > —t±. By inserting (||) into the 
integral (||) one can show that the integral is dominated 
by the constant and by the t± mo+1 terms in @. Car- 
rying out the integral yields Ank F oc t± R ^ \ where the 
exponent i?[mo] approaches quite rapidly its mo — > oo 
limit R = ctjiB — a). For R > 1 this term will again be 
shaded by the linear term as in the case a < ct2 V . In Fig. 
||, we show the dominant and some subdominant expo- 
nents for the spinless and spinful cases as a function of a. 
As one can see, the occurrence of the two-particle regime 
(a > oj2p) reduces the dominant exponent in the spinless 
case, and, in particular, it shifts upwards (from a > 1/2 
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FIG. 2. Exponents controlling the behavior as a func- 
tion of t± of Ank F (thick solid line), and of the FS warp 
AkF (thick dashed line) for the iV-chains system. Thin lines 
indicate subdominant exponents for Arik F : a/(l — a) (solid), 
a/(B — a) (dashed), and unity (dotted). The exponents are 
plotted as a function of the correlation exponent a for the 
spinless and for the spinful cases. 



to a > 2/3) the region where Ank F becomes linear in 
t_i_. It is interesting to compare this result with the pre- 
diction of linear-response theory, which is expected, for 
sufficiently small tj_, either to give the correct behavior 
or to diverge thus signaling a sublinear behavior. The 
linear-response result An k = t±(A + B\k\\ — kp\ 2a ~ 1 ) in- 
deed shows a divergence at kn = kp for a < 1/2, i. e. it 
predicts a sublinear behavior in this region. However, our 
calculation shows that in the spinless case Anj- F is sub- 
linear in t± in a larger region, namely, up to a — 2/3. Un- 
fortunately, this effect is more difficult to see in the case 
of electrons with spin, since in this case ct2 P ~ 0.62 > 1/2 
and the effect is "shaded" by the linear behavior occur- 
ring for a > 1/2. Nevertheless, this effect should be 
detectable, e.g., in the first derivative of Ank F with re- 
spect to k or in the FS warp as discussed below. Notice 
that the exponent R has been obtained by cutting the 
series in t± for the self-energy T _1 at a given order. For 
a > ct2 P , a different choice, like, e.g., cutting the se- 
ries for r, may lead to a different (although in this case 
unphysical) result, which, however, preserves the quali- 
tative effects, in particular the decrease of the exponent 
R. It should be mentioned that recently S. Capponi, D. 
Poilblanc and myself fl2"|| , have evaluated the exponent 
of Ank F by exact diagonalization of small ladders sup- 
plemented by a careful finite-size extrapolation carried 
out by means of an appropriate scaling function. These 
numerical results are in very good accordance with the 
exponents for An^p predicted here and clearly show the 
change of behavior between the two regimes a < ct2 P and 
a > ct2 P - 

We now study the behavior of the Fermi surface 
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as 
fen 



a function of tj_. The FS consists of the points 
H _ given by the solution of the equation |l(]] 
Re Q(ku g \k±],iu = + ;k±)~ 1 = 0. As already dis- 
cussed ]|Jg,§|, within the SDA a solution exists and 
gives fe|| — kp oc ^. At higher order, 

Q{k\\,iio = + ;7r) _1 has the same form as (||) with iui 
replaced with ku — kp and with different coefficients a n , 
whereby B is formally replaced with 1 for a < ct2 P - The 
equation for the (k± — ir) pole can be written in terms 
of a scaling function / as /[tj_/(fc|| — kp) 1 " 01 ] = 0. If 
f[x] =0 has a solution at some point x = x* , then the 
behavior of the FS is similar to the one for the SDA, 
i.e. Ak F = p [tt] - F [0] = (t ± /x*)V(i-«) @. This 
quantity measures the FS warp produced by t±. How- 
ever, it is not easy to find out whether the imaginary part 
of Q^ 1 vanishes fast enough at the FS, so that the above 
result does not necessarily imply that the system goes 
into a Fermi liquid state. Moreover, considering higher- 
order terms for f[x], the solution of f[x] = may not 
exist. In fact, several gaps open in the different modes 
of a TV-chain system fll7|] , which could possibly prevent 
the pole equation from having a solution for some N and 
some k± points of the FS. In this case, as a measure of 
the FS warp, one could define a "pseudo" FS by the point 
fe|| for which \G~ 1 \ has a minimum. Also in this case, the 
FS warp turns out to behave like i^ 1 ^ 1 ""- 1 . Finally, for 
a > ot2p the behavior for small t± is again dominated by 
the dot diagram (a) and by the higher-order diagrams of 
type 7 and 6, as for the calculation of Arik F ■ In a similar 
way, one obtains for the FS warp AkF <x t± R ' a (cf. Fig. 
||). This means that the Fermi surface is warped even 
for a — > 1 in contrast with the SDA result and with the 
expectation coming form the fact that t± is "irrelevant" 
for a > 1. Notice, however, that particle-hole instabil- 
ities producing a gap at the Fermi surface are likely to 
occur for large a. 

In conclusion, using a diagrammatic representation of 
the expansion in power of t±, we have determined the 
power-law behavior as a function of t±_ of the difference 
in occupation An^ F and of the FS warp in a A-chain 
system, our results being exact at any finite order for the 
self-energy. In the single-particle regime (a < ct2 P ), the 
exponents are correctly given by the low-order approx- 
imations. In the two-particle regime (a > ct2 P ) higher- 
order diagrams give dominant contributions and reduce 
the exponents (thus indicating an increased coherence) 
with respect to the SDA result. In this regime, linear- 
response theory in t± is not reliable even in the small-tx 
limit. The increase of the exponents with a, as seen in 
Fig. |^, shows that transverse coherence is reduced with 
increasing interaction but never completely suppressed 
since Ank F and the Fermi surface warp are finite for any 
finite value of t± also for large a. We expect, however, 
that at some value of a, a gap opens due to particle- hole 
instabilities possibly producing incoherence effects in the 
intrachain transport. For a finite number of chains these 



gaps seems to open even for small a |L7| . It would be 
interesting to establish whether a finite critical value of 
a exists in the N — > oo limit. 
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